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Abstract 

Magnetic fields appear wherever plasma and currents can be found. As such, 
they thread through all scales in Nature. It is natural, therefore, to suppose 
that magnetic fields might have been formed within the high temperature en- 
vironments of the big bang. Such a primordial magnetic field (PMF) would 
be expected to arise from and/or influence a variety of cosmological phe- 
nomena such as inflation, cosmic phase transitions, big bang nucleosynthe- 
sis, the cosmic microwave background (CMB) temperature and polarization 
anisotropies, the cosmic gravity wave background, and the formation of large- 
scale structure. In this review, we summarize the development of theoretical 
models for analyzing the observational consequences of a PMF. We also sum- 
marize the current state of the art in the search for observational evidence of 
a PMF. In particular we review the framework needed to calculate the effects 
of a PMF power spectrum on the CMB and the development of large scale 
structure. We summarize the current constraints on the PMF amplitude B\ 
and the power spectral index ns and discuss prospects for better determining 
these quantities in the near future. 
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1. Introduction 



It is well known that magnetic fields influence many physical processes 
over a broad range of scales in the universe. It is natural, therefore, to sup- 
pose that currents and associated magnetic fields could arise from the flow of 
material in the high temperature environment of the big bang. Such a primor- 
dial magnetic field (PMF) would be expected to manifest itself in the cosmic 
microwave background (CMB) temperature and polarization anisotropics 
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and also in the formation of large-scale structure (LSS) [26|, |27j, |28j. These 
studies have pointed out that the effects of a PMF are not negligible if the 
PMF had a field strength B\ ~ 1.0 nG on ~ 1 Mpc scales at the epoch of 
photon last scattering, z ~ 1100. Therefore, it is important to constrain the 
PMF parameters by cosmological observations. The purpose of this review 
is to summarize progress in the development of the calculational techniques 
to deduce such constraints from existing observational data. First, however, 
we motivate this endeavor with a brief discussion of why one thinks there 
might be a PMF from both theoretical grounds and observations of galaxy 
clusters. 



1.1. Generation models of a PMF 

Many authors have been actively studying the origin of cosmological pri- 
mordial magnetic fields. Possible generation mechanism have recently been 



reviewed in detail by |29|. In brief, such primordial fields are expected to 
have random distribution of orientations and field strengths. If the PMF has 
a nearly scale invariant spectrum like the primary scalar spectrum, one of the 
best models is an origin from vector potentials generated during the inflation 
epoch 30, |3l|, |32] . Several authors have proposed that a PMF with a bluer 



spectrum could have been produced during one of the subsequent cosmolog- 



ical phase transitions [33j, |34|, |35|, |36(. A magnetic field could also have been 



generated during or after the epoch of photon last scattering (^1100) by 



vorticities in the cosmological fluid [371 . l38l . l39f . Since each model for the 
generation of the PMF involves different length scales, the spectral index of 
the PMF power spectrum, n B , will depend upon which is the correct model 
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for the generation of the PMF. Therefore, constraining n-Q phenomenologi- 
cally can lead to constraints on the epoch for the generation of the PMF. 



1.2. Observational evidence for a PMF 



Observations of synchrotron emission [40] and Faraday rotation [41], |42 



43| in galaxy clusters all indicate the presence of magnetic fields on large 



scales. Moreover, magnetic fields with a strength of B ~ 1.0 pG have been 



detected in several galaxy clusters |4l|, |42|, |43| . Simulations of such cluster 



magnetic fields may require [44| the existence of a primordial seed magnetic 
field. 

Since the diffusion time of a magnetic field in galaxy clusters is much 
longer than the age of the universe, such a magnetic field is "frozen- in" to 
the ionized fluids [2|. The energy density of the baryon fluid scales as pb oc a~ 3 
while the magnetic energy density scales as pe oc a~ 4 . Therefore, one can 
relate the strength of the magnetic field energy density to the energy density 
of the cosmic baryon fluid -B 3 oc p\. If clusters of galaxies collapse nearly 
isotropically relative to the background space, an observed magnetic field of 
B rsj 1.0 /iG in galaxy clusters now corresponds to a PMF of order ~ 1.0 nG 
at the epoch of photon last scattering near z ~ 1100. 

1.3. Constraint on a PMF by Cosmological observations 

As noted above there is a large difference between the strength of any cos- 
mological magnetic field during the generation epoch and that which could 
be detected at the present time. In most currently proposed theoretical 
scenarios, the generation of a PMF involves some physical process that gen- 
erates current in the early universe before the epoch of photon last scattering. 
However, cosmological observation determine the CMB anisotropies and the 
matter power spectrum on various cosmological scales during and after the 
epoch of photon last scattering. It is therefore timely and desirable to con- 
struct precise theoretical models with which to analyze the CMB anisotropies 
and the large-scale-structure (LSS) matter power spectrum as obtained in 
present and soon to be obtained cosmological observations. The purpose of 
this review is to summarize the development of such models. 



2. Models for the Generation of a PMF 

The origin of a cosmological primordial magnetic field has been an area of 
active research by many authors [See [29| for an excellent review]. The pro- 
posed models can be divided into three classes which we briefly summarize. 
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First, however, note that the cosmological magnetic field damps as 

B oc a~ 2 , (1) 

where a is the cosmological scale factor. In order to avoid confusion, there- 
fore, it is best refer to a comoving PMF field strength Bo scaled to the 
present-day value on some length scale (taken here to be 1 Mpc), i.e. we 
write 



Bo = B{t) x 



air) 



a (r) 



(2) 



where r is the conformal time defined by the differential dr = adt. 



2. 1 . Inflation 

During the epoch of inflation small quantum perturbations are believed 



to have been extended and enlarged into those of super horizon scales [45 



Since spacetime in normal inflation is conformally flat 45], the spacetime 
and any associated electromagnetic field are independent and invariant to 
a conformal transformation. Therefore, unless this invariance to a confor- 
mal transformation is broken, a magnetic field cannot be generated by the 
expan ding spacetime 



30, 31, 32 



In 30 



however, it was assumed that the electromagnetic field and the 
gravitational field interacted in such a way as to break the invariance of the 
electromagnetic fields to a conformal transformation. They then showed that 
a magnetic field could be generated by quantum perturbations. There are 
several means to naturally break such invariance. For example, the strength 
of a PMF generated by an electromagnetic tensor of the dilaton type, has 



been estimated to be ~ InG and ~ 10 nG by [32[ and |46[, respectively. 



2.2. EW transition/ Q CD transition 

Cosmological phase transitions could also produce a PMF with a bluer 



spectrum [33j, |34], |35|, |36| than that of inflation models. For example, bubbles 



of a new lower temperature phase can be made at the cosmological quark- 
hadron 47] or electroweak |48j phase transitions. As those bubbles collide 
and percolate, the baryon symmetry is broken leading to the generation of 



a magnetic field |47J. The strength of the magnetic field generated at these 
epochs has been estimated to be ~ 10~ 7 nG for fields generated by the quark- 
hadron transition 47J and ~ 10~ 14 — 10~ 8 nG for the electroweak transition 

3. 
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2.3. Decoupling and afterward 

Magnetic field generation can also occur on smaller scales during or after 



the epoch of photon last scattering (z < 1000) |38|, |37|, |39[. If there were 
eddies before the recombination epoch, they could generate a PMF 39j. Even 
without turbulent eddies, the known CMB temperature fluctuations imply 
the generation of at least some magnetic primordial field. Since protons have 
mass and photons do not, photons scatter electrons differently than protons. 
This difference induces electric currents as they fall in and out of gravitational 
potentials. These electric fields and can generate a magnetic field of about 



10" 9 nG at 1 Mpc. [39 



3. Constraints on the Physics of a PMF 

In this section we briefly introduce that the effects of a PMF on Big- 
Bang Nucleosynthesis (BBN), the Gravitational Wave Background (GWB), 
the Cosmic Microwave Background (CMB) and Large Scale Structure (LSS). 
We then qualitatively discuss constraints on the PMF parameters deduced 
from each observation. Following that we shall quantitatively review how the 
PMF parameters are constrained by multiple cosmological observations. 

Of course one does not really know a priori the spectrum of cosmological 
magnetic fields. As a starting point, however, it is natural to presume that 
one can characterize the PMF parameters in the same way one does for 
photons, i.e. by a power-law spectrum of PMF energy density fluctuations 
on different wave numbers k, i.e. 

P(B k )<xk nB , (3) 

where P(Bk) denotes the power in magnetic field energy density on the scale 
of wave number k and is a spectral index. 



In previous work {e.g. Refs.p, ll2j} approximate power law spectra were 



utilized to describe a PMF. These approximate spectra have been very con- 
venient for exploring the qualitative effects of a PMF on various physical 
phenomena. However, in this approximation one cannot obtain quantita- 
tively accurate results. For example, one could not accurately constrain the 
parameters of a PMF from cosmological observations. In Yamazaki et al. 
(2008) fl5j. however, a power law formulation for a PMF was developed with- 
out approximation (See Appendix A for details). Using this formula, we 
have been able to calculate the effects of a PMF by accounting for the time 
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evolution of the cut off scale to high accuracy. Also, we have constructed a 
numerical program without approximation, for studying the effects of a PMF 
on the cosmological perturbations. 

Since our main purpose is to review how parameters of the PMF are 
constrained by cosmological observations, we first make a brief explanation 
of the various ways in which a PMF affects cosmology. The reader who 
wishes to understand PMF effects in more detail is referred to the following 
papers regarding the CMB: Ref. @, B H, H Q , LSS: Ref. 0, 0, B H 



and the GWB: Ref. [49j. After reviewing the observational constraints we will 
conclude with a discussion of the constraints upon possible origins of a PMF 
based on the present constraints on the PMF parameters from fits to the 
available observations. 

3.1. Constraints from BBN and the GWB 

The balance between the expansion rate of the universe and various par- 
ticle reaction rates has important effects on the nucleosynthesis of light ele- 
ments in the big-bang. Moreover, since the energy density of the GWB, poW) 
contributes to the total energy density of the universe, the expansion rate 
is affected by the GWB. Therefore, we can indirectly constrain the energy 
density pew from the light element primordial abundances inferred from ob- 



servations of deuterium (D), 3 He, 4 He, and 7 Li. Following 50j one can relate 
the upper limit to the energy density in gravity waves to the upper limit on 
the effective number of neutrino flavors N u present during BBN: 

poo 

/ dtog(i/)/igfi GWB (i/)<5.6xlO- 6 (JV w -3) , (4) 
Jo 

where v is the frequency of the GWB, h is the Hubble parameter (in units 
of 100 km s -1 Mpc^ 1 ), Ogwb is the ratio of the energy density of the GWB 
Pgwb to the critical density of the universe p c . 

The upper limit to the number of neutrino flavors N v is constrained from 
a variety of observations. If one only utilizes a comparison between BBN 
and the inferred primordial abundances of deuterium (D) and helium ( 4 He) 
in the standard big bang one obtains 511 ] : 



N v = S.Utoll at 68% CL 

(2.49 <N V < 3.84 at 68% CL) . (5) 
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Alternatively, if one considers constraints from the WMAP [52|, |53j, |54j, 
CMB power spectrum analysis including Luminous Red Galaxies (LRG) 



55] 
56, 



57] and the Hubble constant, Hq measurements, this limit increases to [58] 

3.45 <N U < 5.01 at 68% CL . (6) 
Or, if one also includes baryon acoustic oscillations (BAO) in the analysis 



this limit becomes 58 



3.46 <N V < 5.20 at 68% CL 



(7) 



An analysis that also includes constraints from Type la supernovae (SN) 



WMAP + BAO + SN + HST leads to 58j 



2.90 <N V < 5.90 at 68% CL 



(8) 



Therefore, depending upon which constraint is adopted, one can rewrite 
Eq.flTJ) as: 



d\og(v)hlQ GWB (v) 



< 
< 
< 



< 4.7 x 10~ 6 BBN, D and 4 He, (9) 

1.13 x 10' 5 (WMAP + LRG + H ), (10) 
1.23 x 10~ s (WMAP + BAO + H ), (11) 
1.62 x 10~ s (WMAP + BAO + SN + HST). 



(12) 



Caprini and Durrer [49] have analyzed the generation of a GWB from 
fluctuations in the PMF. They then used the BBN constraints on the GWB 
to place surprisingly strong constraints on the parameters characterizing a 
PMF. Following 49], we have recalculated the GWB upper limits on the 
PMF field strength as a function of spectral index for a PMF formed during 
different epochs. These are showi£] in FigJTJ The right part of this figure is 
an expanded plot for spectral indices tib from -3 to -2.5. The green, blue, 
and red lines on this plot are the upper limits to a PMF generated during 
big-bang nucleosynthesis, the electroweak transition, or the inflation epoch, 



1 Caprini and Durer (2004) quoted field strength for a comoving scale of galaxies A = O.f 
Mpc, however our results are quoted for galaxy cluster scales A = 1 Mpc. 
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respectively. For this plot we adopt the WMAP constraint on N v from the 



WMAP + BAO + SN + HST analysis [58 . 

If the PMF were generated during inflation, the power should be nearly 
scale invariant (% ~ —3) jij. In this case, the PMF parameters are not con- 
strained strongly by the associated GWB. For example, PMF fields strengths 
as large as B\ ~ 10~ 7 G are allowed for a spectral index of ns ~ —3 (right 
panel of FigJIJ). As we review later, the constraint on the PMF strength from 
various cosmological observations is less than several nano Gauss. Therefore, 
if the PMF were generated during the inflation epoch and the PMF power 
spectrum is nearly scale- invariant, the PMF strength is not constrained very 
much by this method. 

On the other hand, if the PMF were generated after the inflation epoch 
it has a bluer spectrum (smaller spectral index) and is better constrained by 
the GWB. One can deduce from the left panel of FigJIJ that the strength of 
a PMF generated during the EW epoch has an upper limit of order 10 -11 G, 
1(T 15 G or 10~ 20 G for n B = -2, -1 or 0, respectively. Similarly, if the PMF 
were generated during the BBN epoch, the upper limit of the PMF strength 
is constrained to be order 10 -10 G, 10~ 14 G or 10 -18 G for n-Q = -2, -1 or 0, 
respectively. 

Furthermore, since the upper limit on the PMF amplitude, B^ p , depends 
upon the effective number of neutrinos N v allowed in the early universe: 

Bf oc y/W u - 3, (13) 

the constraints on the PMF parameters strongly depend upon independent 
observations which constrain N v . 

There has also been an attempt to constrain the PMF field strength from 
direct measurements of limits to the present cosmological GWB. |59] has 
deduced the constraint B A <4x 10" 7 G at 1 Mpc and B A <9x 10~ n G at 
100 Mpc from the LIGO S5 data . 



3.2. Constraints from the CMB 

For a PMF of order 1.0-10 nG at the surface of photon last scattering, the 
total energy density in the PMF is much smaller than that of the temperature 
fluctuations of the CMB (called the T-mode, where T denotes the total scalar 
temperature). Therefore, we can treat the energy density of a PMF as a first 
order perturbation and assume a stiff source for the time evolution . That 
is, all back reactions from the fluid onto the magnetic field can be discarded 
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because they are second order perturbations. In this case, we can also assume 
that the conductivity of the primordial plasma is very large and that the 
electric field is negligible, i.e. E ~ 0. This "frozen-in" condition is a very 
good and useful approximation 0. 

Recently, the CMB polarization power spectra have been observed by 



several projects. This has allowed for more precise constraints [231 . 1611 . 162 
on the PMF and other cosmological parameters. For clarity, we first give a 
brief description of the CMB polarization as follows: A photon scattered by 
an electron is polarized perpendicular to the incident direction. When the 
incident photons are isotropic or have only a dipole distribution, there is no 
net polarization of the scattered photons. On the other hand, if the incident 
photons are perpendicular to each other and have different intensities, the 
scattered photons will have a net linear polarization. In addition, photons 
are polarized by the perturbed gravitational potential, e.g. weak lensing 
effects. The polarization of photons with positive parity is called the " E- 
mode", and the polarization with negative parity is called the "5-mode". 
Therefore, there are three observable modes: the T-mode; the .E-mode and 
the S-mode. From these, one can generate correlations leading to three power 
spectra denoted: TT, EE, and BB. There are also three cross-correlation 
spectra: TE, TB, and EB. The only nonzero spectra, however, are the TT, 
EE, BB, and TE modes due to parity considerations. Also, the modes can 
have up to three kinds of fluctuations: scalar, vector and tensor. 

Figure [2] illustrates that a PMF produces the largest temperature fluctua- 
tions and polarization anisotropics of the CMB for higher multipole moments 
i. The energy density of the PMF is proportional to the fourth power of the 
scale factor a -4 just like the photon energy density. On the other hand, the 
field strength of the PMF scales with the plasma density p p0 as 

B\(k, a) = B x {k) * {p p0 + 5p p ) (14) 

where 5p p is a first order perturbation in the plasma density. If the field 
strength of the PMF is less than the order of 10 nG, the energy density of 
the PMF is less than or equal to the perturbations in the energy density 
of the CMB photons. Thus, the second term of Eq. (jl4p is a second order 
perturbation and negligible in the linearized theory. Therefore, we can ignore 
the effect of the time evolution of the plasma density fluctuations (e.g. Silk 



damping |63|) on the PMF in the linear theory. As a result, the effects of 
a PMF on the temperature fluctuations and polarization anisotropies of the 
CMB tend to be large even for scales smaller than the Silk damping scale. 
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For the vector mode, the CMB polarization due to a PMF is largest for 
higher multipole moments, £, while the tensor mode of the PMF diminishes 
for large multipoles [cf. Figj2] [Hi], 0, 15, 23]. The reason for this is that the 
gravity waves from the PMF can be negligibly small after horizon crossing. 
They are small because, once inside the horizon, the homogeneous solution 



for the gravity waves begins to oscillate and decay rapidly. [65|, |66|, |67|, |68[ . As 
a result, the effect of gravity waves from the PMF only occurs on scales larger 
than the horizon at the epoch of the generation of the CMB (z ~ 1000). 

There are several degeneracies between the PMF and other physical pro- 
cesses in the early universe. Panel (b) on Figures |2] and |3] shows the BB mode 
of the CMB polarization from the primary fluctuations, gravitational lensing, 
and the PMF. We can see that the total BB mode spectrum is dominated 
by the PMF for B\ > 2.0 nG and £ > 200. We must, however, consider 
the degeneracy between the PMF and lensing effects on such small angular 
scales. This is because the BB mode is converted from the EE mode by a 



gravitational lensing effect on these angular scales [681 ] . However, the source 
of the spectrum from the gravitational lensing signal is the EE-mode of the 
CMB polarization of the primary scalar fluctuations. Fortunately, however, 
this spectrum can be subtracted directly because the E'.E'-mode polarization 
of the primary scalar fluctuations have been determined independently. 

Since the Sunyaev-Zel'dovich effect is one of the main foreground sources 
for higher £ [69], we must also consider the degeneracy between the PMF and 
the Sunyaev-Zel'dovich effect on the temperature fluctuations of the CMB. 
Nevertheless, the effects of a PMF on the CMB background are independent 
of frequency and the Sunyaev-Zel'dovich effect on the CMB as a foreground 
depends upon the frequency. Therefore, it should be possible to distinguish 
the effects of a PMF from such foreground effects by making observations at 
different frequencies. 

The qualitative features of the constrained parameters of the PMF from 
only the CMB observations are easy to understand. Table [1] shows the con- 
straint on the PMF parameters from the combined CMB observations of 
WMAP, ACBAR and CBI using Markov Chain Monte Carlo (MCMC) meth- 
ods [70]. We find that there is no obvious degeneracy between the PMF and 
cosmological parameters. The reason for this is that the effects of a PMF 
dominate for £ > 1000, while the other cosmological parameters are con- 
strained by the WMAP power spectrum obtained for £ < 1000. 

On the other hand, from Fig. H]one can see that there is a strong degen- 
eracy between the PMF amplitude B\ and the power law spectral index n-Q. 
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There are two main reasons for this. For one, the effect of a PMF on the 
CMB is to produce a peak in the range of £ > 1500. Unfortunately, however, 
observations are not yet precise enough in this multipole range to constrain 
the PMF parameters. This is especially true for the BB mode for higher 
I where we do not yet have observations of the BB mode with sufficiently 
small error. The second reason is simply that the PMF parameters them- 
selves induce similar effects on the amplitude of the CMB power spectrum 



12, 15 



3. 3. Constraints from LSS 

There are strong constraints on parameters of the PMF from the observed 



limits on the erg parameter [see 28J for details] . This parameter is determined 



from a weighted integral of the matter power spectrum [71| and corresponds 
to the root-mean-square of the matter density fluctuations in a comoving 
sphere of radius 8/i _1 Mpc. 

Figure. [5] shows that the density fluctuations of matter are more strongly 
affected by a PMF for wavenumbers k/h > 0.1 Mpc _1 [28|]. As mentioned in 
Sec. I3.2[ the time evolution of the PMF energy density does not depend on 
the time evolution of the plasma density fluctuations (e.g. Silk damping) in 
the linear theory. Therefore, the PMF can survive as a source of temperature 
fluctuations and polarization anisotropies of the CMB on scales well below 
the photon diffusion length. 

The baryons influence the cold dark matter(CDM) through gravity. This 
effect is very small before the epoch of photon last scattering (z ~ 1000) 
because the baryon density oscillates with the photons and their gravitational 
effect on the CDM becomes very small. However, after the baryons decouple 
from the photons, the baryon density begins to affect the density of the CDM 
through gravitational interactions [72j. Also, the baryons are influenced by 
any PMF present in the early universe. Thus, a PMF can indirectly affect 
the CDM evolution. Nevertheless, it is reasonable to assume that the energy 
density of the PMF is independent of the matter density fluctuations. In this 
case, the PMF merely increases the matter power spectrum independently 



of whether the pressure or tension dominate the PMF [28 



As discussed in Sec. 13. 2^ there is the strong degeneracy between the PMF 
amplitude B\ and the power law spectral index n-Q. Thus, one needs a good 
Bayesian prior constraint to effectively decouple these two parameters char- 
acterizing the PMF. In this regards one is aided by the observed constraints 
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on the erg parameter. The parameter <r 8 is constrained by observational data 
on linear cosmological scales [73|, 12^, |57|, |56| to lie in the range 0.7 < cr 8 < 0.9. 

Fortunately, the recent CMB observations have determined cosmologi- 
cal parameters based upon the power spectrum obtained on larger scales 

The effects of the PMF are mainly on smaller scales ( 
Therefore, one expects that there is only a small degen- 



{£ < 1000H54]5l]53 



12 



28. 15 



£ > 1000) 

eracy between the PMF parameters and the other cosmological parameters. 
Hence, if one wishes to understand the qualitative nature of how the PMF 
parameters are constrained by observations of LSS, one is justified in fixing 
the other cosmological parameters at their best fit values from the WMAP 
analysis. 

In Fig. O we show the optimum PMF parameters n-Q and B\ for various 



constant values of as- The power spectrum of the PMF PpMF{k) scales 15 
as: 



PMF 



(k) oc k 



2n B +3 



(15) 



Therefore, for n-Q < —1.5, the effect of a PMF on the density fluctuations 
on small scales decreases with smaller values for n-Q. In the case that ub is 
-3.0, the matter power spectrum including the PMF effects diminishes 



near 



for smaller scales. Hence, larger amplitudes of B\ are allowed. However, for 
tib > —1.5 the PMF power spectrum is proportional to the cut-off scale as 



t 2n B +3 



15 ] where kc is proportional to B x 1 




Substi- 



tuting these relations into Eqs.(15)-(17) in Ref. [15[, an expression for the 
relation between the power spectrum of the PMF and the magnetic strength 
can be deduced, 



P(k) 



PMF 



oc B 



[14/(n B +5)] 



(16) 



When B\ « InG, the matter power spectrum including the effect of the 
PMF for n-Q > —1.5 becomes larger for larger tib- In this case a large strength 
of B\ is not allowed for larger ub- 

If the PMF has no correlation with the primary density fluctuations, the 
effect of a PMF is to increase the matter power spectrum independently of 



whether the PMF pressure or tension dominates [15] . As noted above, recent 



cosmological observations of LSS imply erg > 1 can be excluded. Therefore, 
we can exclude field strengths of B\ >1 nG if ub > —0.9 (Fig. [6]) and B\ >0.1 
nG if ub > 0.2. When the PMF evolves to the observed magnetic field in 
clusters of galaxies, we can expect that the PMF amplitude is Bx j^l nG. In 
this case from FigJHlwe can exclude spectral indices of tib > —0.9. 
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3.4- Concordance MCMC analysis 

Constraints on parameters of the PMF have been deduced by our group 



231 ] and other groups [61|, |62|- These constraints are based upon fits to the 
CMB and LSS observational data in the context of a flat Lambda CDM 
cosmology. This cosmology is characterized by six standard parameters 
i.e. {Vt b h 2 , Vt c h 2 , t c , n s , log(10 10 v4 s ), A t /A s , \B\\, n B }, where Vt b h 2 and Q c h 2 
are the baryon and CDM densities in units of the critical density, h de- 
notes the present Hubble parameter in units of 100 km s Mpc , tc is the 
optical depth for Compton scattering, n s is the spectral index of the primor- 
dial scalar fluctuations, A s is the scalar amplitude of the primordial scalar 
fluctuations and A t is the scalar amplitude of the primordial tensor fluctu- 
ations. We define the tensor index for the primordial tensor fluctuations as 
n t = — (A s /A t )/S. For these standard parameters one can adopt the Bayesian 
priors used in the previous WMAP analysis 54|, |52|, |53|, |78[ without a PMF. 
To the standard cosmological parameters, one then adds the two PMF pa- 
rameters: the spectral index n B ; and the average magnetic field strength on 
a comoving scale B\. 

Using a MCMC method with cosmological observations (e.g. the CMB 
and/or LSS) the standard cosmological parameters and the PMF parameters 
have been constrained [23J. The results of this analysis are summarized in 
Table [TJ However, we note that an inherent flaw in the MCMC approach is 
a sparse sampling of points near the boundary. Consequently, this method 
tends to assign low probability near a boundary. This is the reason that 
the PMF power law index n# seems to be constrained even in the limit of 
B\ — > 0. In reality, of course, the spectral index ns is not constrained for 
B\ = 0. However, this shortcoming does not negate the result that a finite 
magnetic field and spectral index give a genuine minimum likelihood in the 
goodness of fit. In the best fit including the PMF parameters, the minimum 
total x 2 is improved from 2803.4 to 2800.2 corresponding to an improvement 
of the reduced \ 2 from 1.033 to 1.031. This slight improvement implies that 
the existence of the PMF is still only of marginal significance. 

Figures [7] and E] show the 68% and 95% C.L. probability contours in the 
planes of each standard cosmological parameter versus field strength of the 
PMF or power law index. Also shown are the probability distributions. The 
bottom panels of Fig. [7J and Fig. [8] show the probability distributions for 
\B\\ and n B - Of particular note for this MCMC analysis is that a parameter 
set of \B\\ = 0.85 ± 1.25 nG and n B = —2.37^'^ maximizes the likelihood. 
These values of the PMF parameters are consistent with no magnetic field, 
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and thus provide only imply upper limits. Nevertheless, they suggest that it 
may be possible to detect a finite PMF with forthcoming data (particularly 
for large CMB multipoles). 

These figures also confirm that there is no degeneracy between the stan- 
dard cosmological parameters and the PMF. We can see this from Table [1] 
in which it is apparent that the standard cosmological parameters do not 
significantly differ from those deduced directly from the WMAP data with- 
out the introduction of a PMF. As mentioned above, the reason for this is 
simple. The standard cosmological parameters are mainly derived from the 
observed CMB for low multipoles t < 1000 (up to the 2nd acoustic peak), 
while the PMF is most important for i > 1000. Thus, the PMF effect on 
the CMB power spectrum is nearly independent of the standard cosmological 
parameters. 

The tensor to scalar ratio A t /A s in our analysis is defined as the primary 
tensor amplitude without the PMF. When we compare our tensor amplitude 
with the cosmological observations by the MCMC method, we combine the 
tensor amplitude from the PMF A^pmf] and A t . Therefore the upper limit 
to the tensor to scalar ratio A t /A s in our analysis is less than the A t /A s 
ratio constrained from the WMAP analysis without a PMF (A t /A s < 0.43, 
95% CL). Actually, the upper limit of (A t + A t [ PMF })/A s constrained by our 
analysis including a PMF is consistent with the previous constraints. In 
addition, the effects of a PMF on the matter power spectrum and the CMB 
temperature fluctuations and polarization anisotropies dominate for smaller 
scale e.g. i > 1000, while the PMF contribution to the CMB for larger scales 
is negligible compared to the primary CMB fluctuations of the scalar and 
tensor modes (l6l . 79]. Therefore, the tensor to scalar ratio is not affected by 
the presence of a PMF. 

Figure [S] shows the deduced probability distributions and the la and 2a 
(68% and 95% C.L.) probability contours for the derived parameters, as, H , 
-Zreion, and Age, where H is the Hubble parameter in units of km s _1 Mpc -1 , 
z r eion is the red shift at which re-ionization occurs, and Age is the presently 
observed age of the universe in Gyr. In Fig. [5J one can see by the slight bend 
in the contours that there is a weak degeneracy between the PMF parameters 
and as- The reason for this is that the matter power spectrum for smaller 
scale determines both as and the PMF parameters. 

Table [1] summarizes the constrained PMF parameters together with the 
standard and derived cosmological parameters. The best constraints on the 
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PMF determined yet to date determined in our analyses are [2i 

|B A |< 2.10 nG (68%CL) , < 2.98 nG (95%CL) (17) 
on a present scale of 1 Mpc, and 

n B < -1.19 (68%CL) , < -0.25 (95%CL). (18) 



This result differs only slightly from the constraints obtained in [80|, |62J based 
upon different observations. 

Adopting these PMF constraints and the BBN GWB constraint discussed 
above ill, a re-examination of the three main PMF generation scenarios 



33l . l34j . l35l . |36| is motivated by the present results. Figure [9] summarizes the 
contours of \B\\ and hb at the 1 and 2 a CL. along with the associated 
GWB constraints for the three generation epochs. 

The region bounded by the yellow area corresponds to the 2a CL. on 
PMF parameters. The upper limit of the produced PMF from big-bang 
nucleosynthesis, the electroweak transition, and the inflation epoch are shown 
by the green, blue, red lines, respectively. However, as mentioned above, 
there is not yet a lower limit for the PMF parameters, B\ and n B at the 
2cr C.L. On the other hand, considering the region bounded by the a C. 
L. red orange area, one can find both upper and lower limits to the PMF 
strength. Figure. [9] shows that the allowed or excluded regions based upon 
these multiple constraints depends upon when the PMF was generated as 
followsjfl 

At the 2a C.L. there are only upper limits to \B\\ and no limits on 
These are: 



I. Inflation: 

\Bx\ < 2.90 nG 



(19) 



II. Electroweak transition: 

\B X \ < 3.08 nG 



(20) 



2 The 1 a and 2a regions on the two parameter plane like Fig|9]do not directly indicate 
constrained values of the parameters at la and 2a as given in Tabled] 
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III. BBN: 



\B X \ < 3.10 nG 

(21) 

At the la CL. there are both upper lower limits to \B X \ and n B . These 

are: 

I. Inflation: 

0.292 nG < \B X \ < 2.33 nG 
-2-97 < n B < -2.66 

II. Electroweak transition: 

0.117 nG < \B X \ < 2.46 nG 

-2-97 < n B < -2.11 

III. BBN: 

0.117 nG < \B X \ < 2.48 nG 

-2-97 < n B < -2.03 

If the PMF were generated at an even earlier epoch, it is clear that upper 
limits on both \B X \ and n B would be more stringent. These limits are the 
strongest constraints on the PMF that have yet been determined. However, 
one should keep in mind that the evolution and/or the generation of the 
magnetic field on cosmo logical scales during the formation of LSS is not well 
understood. Therefore, if there are other effective physical scenarios for the 
generation and/or evolution of the PMF after the epoch of the photon last 
scattering, our lower limits to the PMF parameters will change. To accurately 
constrain the PMF parameters, one should study the PMF not only before 
but also after the epoch of photon last scattering. 

4. Discussion and Future Challenges 

The generation and/or the evolution of the PMF along with its physical 
behavior in the early universe have been studied by many researchers [see 



29j for a recent review]. In this section we review several relevant topics 
and their implications for different PMF strengths. We also discuss future 
prospects for studying the role of the PMF in cosmology, astrophysics and 
astronomy. 
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Cosmological Parameters 



Parameter 


mean 


best nt 


Q b h 2 


0.02320 ± 0.00059 


0.02295 


Q c h 2 


0.1094 ± 0.0046 


0.1093 


rc 


U.Uo ^ ± U.U1 ( 


n no o 


n s 


0.977 ±0.016 


0.970 


ln(10 10 A s ) 


3.07 ±0.036 


3.06 


A t /A s 


< 0.170(68%CL), < 0.271(95%CL) 


0.0088 


\B x \(nG) 


< 2.10(68%CL),< 2.98(95%CL) 


0.85 


n B 


< -1.19(68%CL),< -0.25(95%CL) 


-2.37 


0"s 


n oio+U.U28 
U.O±Z_ 033 


0.794 


H 


73.3 ±2.2 


72.8 


^reion 


10.9 ± 1.4 


10.5 


Age(Gyr) 


13.57 ±0.12 


13.62 



Table 1: PMF parameters and ACDM model parameters and 68% confidence intervals 
(A t /A s is a 95% CL) from a fit to the seven year WMAP 81] + ACBAR [H + CBI [H 
+Boomerang 84] + 2dFDR [z| data. 

4.1. The BB mode from the PMF 

As mentioned above, a PMF affects not only the temperature fluctuations, 
but also the polarization anisotropics of the CMB. Although fits have been 
based upon all available temperature and polarization data of the CMB, it 
turns out that the TT and BB modes are the most important for constraining 
the parameters of the PMF. Figs. [10] (the TT-mode) and [TT] (the BB-mode) 
show a comparison of the observations of the CMB with the computed best- 
fit total power spectrum. Figure [10] shows the best fit and allowed regions 
including both the Sunyaev Zel'dovich (SZ) effect (scattering from re-ionized 
electrons) at the K(22.8GHz) band (upper curves) and without the SZ effect 
(lower curves). Including the SZ effect only slightly diminishes the best fit 
magnitude of the PMF. 

Figure HU] exhibits the best-fit and allowed regions (or curve) both in- 
cluding the SZ effect at 22.8GHz (upper curves) and without the SZ effect 
(lower curves). Since the temperature fluctuations of the CMB from recent 
observations have become more accurate than the data shown in Fig. [10] 
these observations give a strong prior range for constraining parameters of 
the PMF. However, as discussed above, we cannot determine the parameters 
of the PMF by the TT mode alone because there is a strong degeneracy 
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between the field strength and the power law index of the PMF. On the 
other hand, the effects of a PMF on the BB-mode of the CMB has a feature 
for i > 200 as shown in Fig. [TTJ This is due to the vector source from 



the PMF|12|. |15| . Thus, if accurate observations of the BB mode could be 
obtained in the future, one expects that the parameters of the PMF could 
be constrained without the uncertainty due degeneracy between the field 
strength and the power law index. 

4.2. PMF and GWB 



Several groups |49|, |59] have studied the GWB induced by a PMF. From 
such studies, one expects that there is a degeneracy between the GWB from 
inflation and that generated by a PMF. Indeed, if the PMF is generated 
during the inflation epoch this degeneracy would become stronger. This is 
particularly true since the power spectrum of an inflation-generated PMF is 
expected to be nearly scale- invariant. Therefore, even if the GWB is observed 
in the future, it will still be difficult to constrain the parameters of the PMF 
from these observations alone(See Fig. 12). 

There is, however, a difference. The GWB from inflation mainly affects 
the tensor BB mode of the CMB, while a PMF affects both the vector and 



tensor polarizations of the BB mode of the CMB [15|, |23j. Thus, the BB 
mode exhibits different features in the GWB from a PMF or an inflation 
origin. Since both a PMF and inflation dominate the tensor source on the 
same scales of the BB mode, there is a strong degeneracy between the BB 
mode of tensor type polarizations from the PMF and inflation. However, 
the PMF also has a vector type polarization as a source of the BB mode 
of the CMB and this mode is most important on smaller angular scale e.g. 
I > 200, while both tensor sources of the PMF and inflation dominate the 
BB mode on larger angular scales e.g. i > 200. Therefore, if one could 
simultaneously obtain observations of the GWB and the tensor and vector 
BB modes of the CMB, this degeneracy may be resolved. 

4-3. Constraint on the neutrino mass in the presence of a PMF 

The presence of a PMF can alter the particle constraints deduced from the 
CMB. If the velocities of finite-mass neutrinos are large enough, the growth 
of density fluctuations is impeded on the free-streaming-scale of the neutrinos 



85| . Hence, the matter density field in the early universe can be affected by a 



finite neutrino mass |86|. The inferred cosmological constraints on the mass 



of the neutrino are of considerable interest. The sum of the neutrino masses 
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is constrained to be V) m v < 0.1 — 1 eV by a combination of cosmological 
observations 1871 . l88l . l89l . l58| and the results of tritium beta-decay endpoint 



experiments [90(. On the other hand, as explained Sec. I3.3[ a PMF affects 
the matter density fields in the early universe and the parameters of the 
PMF correlate with o~%. Therefore, there is a possible degeneracy between 
the parameters of the PMF and the mass of the neutrino. This degeneracy 
relaxes the cosmological constraint on the sum of the neutrino masses. 

In Fig. [13] we show the constraints on the sum of the neutrino masses 
~Y^v m v an d B\ for various fixed values of tib and ranges of erg. The expected 
field strength of the PMF from cosmological observations is B\ <2.0nG(lcr) 
and <3.0nG(2er), while erg is constrained to be in the range 0.75< o~$ <0.85. 
For this range of <7g, if we do not consider the effect of a PMF on the matter 
density field, the constrained sum of the neutrino masses is less than 0.11 
eV. On the other hand when a PMF is introduced, the upper limit on the 
sum of the neutrino masses from Fig. [131 increases to m v < 0.24 eV 
for Ub = —1.5 and < 0.6eV for = —2.5 and N v = 3. Since the effect 
of a PMF cancels the effect of neutrinos on the density fluctuations, the 
constrained mass of the neutrinos when a PMF is present is larger than that 
deduced without considering a PMF. 

However, <jg depends upon the cosmological model employed and it has 
some degeneracy with fl m = Q b + fi CDM , an d A s , even if these pri- 



mary parameters are well constrained by the CMB data, e.g^ WMAP 58 



Furthermore, the neutrino mass has a degeneracy with Vt m 91]. Thus, one 
should ultimately consider the simultaneous degeneracy between the stan- 
dard cosmological parameters, the mass of neutrinos and the PMF. With 
high precision observations on small angular scales, it may be possible in the 
near future to obtain not only upper but also lower limits to the mass of the 
neutrinos from cosmology and astrophysics in the presence of the PMF. 

4-4- Non-Gaussianity from a PMF 

If we assume that the inflation mode is simple, then the primordial fluctu- 
ations will be Gaussian. A measurement of non-Gaussianity in the primordial 
fluctuations could thus provide important evidence for new physics beyond 
that of simple inflation [92]. Recently, however, several authors have shown 



that a PMF can cause non-Gaussianity [93|, |94], [95|, |96|, |97 . 

We can constrain the non-Gaussianity by auto- and cross correlated bis- 
pectra from the intensity fluctuations of the CMB. On the other hand, the 
bispectra from the CMB are also affected by a PMF. Therefore, one expects 
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that the effects of a PMF and inflation on the non-Gaussianity of the CMB 
will have a degeneracy. Hence, even if one could obtain accurate observations 
of the non-Gaussianity of the CMB, one would not be able to analyze the 
inflation model without information on the effects from the PMF. To solve 
this degeneracy and truly understand the inflation model, one should con- 
sider the constraints on PMF effects by various observations, e.g. BBN, the 
GWB, LSS, and the CMB. 

4-5. 21 cm line and a PMF 

The 21 cm line is emitted by a transition of the neutral hydrogen atom 
between the two different energy levels of its Is ground state. Since hydrogen 
is the most abundant element in the Universe (75 % of normal matter by mass 
and 90 % by number), it is important to survey the 21 cm line for researching 
the "Dark Age'fl If observations of the 21 cm are successful, one should be 
able to obtain a precise matter power spectrum after recombination and the 
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As mentioned above, a PMF affects structure formation in the early Uni- 
verse. Furthermore, a PMF of strength ~ 1 nG reheats the intergalactic 
medium (IGM) including hydrogen up to several thousand K through am- 
bipolar diffusion [19j. Therefore, the 21 cm line is indirectly affected by a 
PMF. The Square Kilometer Array (SKA) @ is planning a most promising 
future project for observing the cosmological 21 cm line. One should con- 
sider the effects of a PMF on the 21 cm line to obtain an accurate physical 
interpretation from the results of observations like that of the SKA. 



3 Considering the cosmological redshift, the 21 line is observed in a frequency range 
from about 200MHz to 9 MHz. 

4 the detected frequency range of the SKA project will be from 70 MHz to 10GHz. 
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Inflation EW BBN 

WMAP+BAO+SN+HST(Upper limit 68%CL) 

the BBN theory(Upper limit 68%CL) 

the BBN theory(best fit) 




Figure 1: Constrained parameter plane of the PMF field strength B\ vs. the power 
spectral index n-Q by the BBN. The thin (green), middle (blue), and bold (red) lines show 
upper limits on the produced PMF during the epoch of BBN, the electroweak transition, 
and the inflation epoch, respectively. (For interpretation of the references to colour in this 
figure legend, the reader is referred to the web version of this article.) 
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Figure 7: Probability contour plane of PMF parameters (n-g and B\) vs. the standard 
cosmological parameters and probability distributions for PMF parameters and the stan- 
dard cosmological parameters Deep (red) and pale (orange) contours show the 1 a(68%) 
and 2 cr(95%) confidence limits. Bold (red) curves on the bottom of the figure are the 
probability distributions for the PMF parameters. Dashed (black) curves on the right of 
the figure are the probability distributions of each standard cosmological parameter. (For 
interpretation of the references to colour in this figure legend, the reader is referred to the 
web version of this article.) 

35 



\B\\ [nG] n B Probability 



1 2 3 -3 -2 -1 0.2 0.4 0.6 0.8 1 




1 2 3 -3 -2 -1 0.2 0.4 0.6 0.8 1 

1^1 [nG] n B Probability 



Figure 8: Probability contour plane of PMF parameters (ns and B\) vs. the derived pa- 
rameters 08, Hq, z r eion and Age and probability distributions for the derived parameters. 
Deep (red) and pale (orange) contours show the 1 <r(68%) and 2 cr(95%) confidence limits. 
Dashed (black) curves on the right of the figure are the probability distributions for each 
derived parameter. (For interpretation of the references to colour in this figure legend, the 
reader is referred to the web version of this article.) 
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Figure 9: Allowed and excluded regions in the la (68%) C.L. and 2a (95.4%) C.L. in the 
plane of IBi | vs. n B obtained by an MCMC method applied to the CMB (WMAP 5yr 0, 
ACBAR[i3], CBl[83|, Boomerang 0]) and the LSS (2dFDR j73|) observational data. The 
thin (green), middle (blue), and bold (red) lines show upper limits on the produced PMF 
during the epoch of BBN, the electroweak transition, and the inflation epoch, respectively. 
(For interpretation of the references to colour in this figure legend, the reader is referred 
to the web version of this article.) 
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Figure 10: Comparison of the observed CMB with the computed total TT power spec- 
tra from the best-fit and allowed parameters. Various ranges of each plats are (a) 
TT(2 < £ < 2750) and (b) TT(400 < I < 2750). Deep colour (red) regions indicate 
the best-fit parameter set and allowed regions [pale (orange) colour] are from the con- 
strained parameter set of Table. 1. Dots with error bars show the WMAP 5yr, ACBAR 
08, Boomerang and CBI data as the legend box in the each panel. The upper curves of 
allowed regions in each panel include the SZ effect at the K(22.8GHz) band, and the lower 
curves do not include the SZ effect. (For interpretation of the references to colour in this 
figure legend, the reader is referred to the web version of this article.) (For interpretation 
of the references to colour in this figure legend, the reader is referred to the web version 
of this article.) 
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Figure 11: Comparison of the observed CMB with the computed total BB power spectra 
from the best-fit and allowed parameters. Deep colour (red) regions indicate the best-fit 
parameter set and allowed regions [pale (orange) colour] are from the constrained param- 
eter set of Table. 1. Dots with error bars show the WMAP 5yr, ACBAR 08, Boomerang, 
CBI, CAPMAP, and DASI data as the legend box in the figure. Downward arrows for 
the error bars indicate that the data points are upper limits. (For interpretation of the 
references to colour in this figure legend, the reader is referred to the web version of this 
article.) 
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Figure 12: Matter power spectra from finite mass neutrinos and a PMF. The range of 
the left and right panel is 0.01Mpc~ 1 k/h < 0.2Mpc _1 and 0.1Mpc _1 fc/ft, < 0.2Mpc _1 
Dots with error bars and curves are as indicated in the legend box in the figure. To better 
illustrate effects of finite-mass neutrinos and a PMF for the matter power spectra, we have 
used a larger field strength of a PMF B\. This figure shows that a PMF affects strongly 
the matter power spectra on k/h > O.lMpc -1 . The magnitude of PMF effect depends 
on both B\ and n^. On one hand, the total amplitudes of the matter power spectra are 
decreased by the neutrino mass. Furthermore, the effects from the mass of neutrinos on 
larger k (smaller scales) are greater than on smaller k (larger scales). 
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Figure 13: Allowed and excluded regions for the as on the B\ vs. J2n =3 m v plane. Filled 
regions indicate ranges of a% as 0.75 < <r$ < 0.85 and pale colour (gray), middle colour 
(yellow) and deep colour (red) regions show B\ > 3.0 nG, B\ < 3.0 nG and B\ < 2.0nG, 
respectively. (For interpretation of the references to colour in this figure legend, the reader 
is referred to the web version of this article.) 
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